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Abstract 

We construct the theta function solution to the initial value problem 
for the generalized periodic discrete Toda equation. 

1 Introduction 

The aim of the present paper is to obtain an explicit formula for the solution to 
the hungry periodic discrete Toda equation (hpdToda) (|1.1H1.3|) : Vn., t £ Z, 

i n +M = i n + v*-v*±{, (i.i) 

rt yt 
V n - T t+M ' 

I n ~ In+N> Vn — Ki+JVj (1-3) 



where N and M are positive integers, t is the time variable and n means the 
position, and relation (|1.3|) is just the periodic boundary condition. This system 
is a variant of the periodic discrete Toda equation, which is the M = 1 case [4]. 



This article is a continuation of the paper pp. We will construct a tau 
function solution for the hungry periodic discrete Toda equation (hpdToda). 

Remark: To avoid a non- interesting solution I n +M = V^ +1 = In+i, we 
should assume the extra constraint 

y\ N Tt+M _y\ N 7* V\ N V t+1 - Y\ N V 1 

lln=l J ra — lln=l 1 n 7~ lln=l v n ~ lln=l V TI1 

which is enough to guarantee the existence of a unique solution. See theorem 
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Notation: For a meromorphic function / over a complete curve C, (/)o 
(resp. (/)oo) denotes the divisor of zeros (resp. poles) of /. Let (/) := (/) — 
(/)oo- Div d (C) means the set of divisors over C of degree d and Pic d (C) 
means the quotient set defined by Pic d (C) = Div d (C) /(linearly equivalent). 
For an element V G Div d (C), [V] means the image of T> under the natural map 
Div rf (C) ^Pic d (C). 

In sections [2] and O we consider the case g.c.d. (N,M) = 1. Wc will discuss 
the general cases in section SJ 



2 Linearization of hpdToda 

We summarize the results of [I] briefly in this section. The reader should consult 
the paper for further details. 

2.1 The spectral curve and the eigenvector mapping 

The hpdToda equation (|1.1H1.3|) is equivalent to the following matrix equation: 



L t+1 (y)R t+M (y) = R t (y)L t {y), 
where L t {y) and Rt{y) are given by 

. , Mv) = 
i 



(2.1) 



Lt(y) 



( l 



i 



v 



v JV-l 



/ 



P 



\y 



i 



and y is a complex variable. Let us introduce a new matrix X t (y) defined by 

X t (y) :=L t (y)R t+ M-i(y)---Rt+i(y)Rt(y). (2.2) 

From (j2~Tj) and (|2~2|) . we obtain 

X t+1 (y)R t (y) = R t (y)X t (y), (2.3) 

which implies that the characteristic polynomial of X t (y) is invariant under 
the time evolution. Let F(x,y) :— det {X t (y) — xE) be the characteristic poly- 
nomial of X t (y) (E is the unit matrix). Denote the afiine curve defined by 
F(x, y) = by C, and its completion by C. Of course, C is invariant as well 
under the time evolution. This projective curve C is called the spectral curve 
of the hpdToda. 
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2.1.1 Properties of the spectral curve 

Now let us list the behaviour of C, following [T] §2. 

• on C, there exists a point P : (x, y) — (00,00) around which there exists 
a local coordinate k such that x = k~ M + ■ ■ ■ and y = k~ N + ■ ■ ■ . 

• on C, there exists a point Q : [x, y) = (00, 0) around which there exists 
a local coordinate k such that x — Ek^ 1 + ■ ■ ■ and y — k N + • • • , where 

• the M points Aj : (x, y) = (0, (-1)^ \[ n P n ), j = 0, 1, . . . , M - 1 lie on 
C. 

• the point B : (x,y) = (0, lies on c - 

• The projection p x : C 9 (x, y) 1— > a; G P 1 is (M + 1) : 1, and the projection 
p y : C 9 (ar, y) 1-* y e P 1 is N : 1. 

• C has genus g = ( jv ~ 1 )( M + 1 )~ m + 1 ^ w here m is the greatest common divisor 
of N and M. 

Hereafter we assume C is smooth unless otherwise stated. 

2.1.2 The eigenvector mapping 

An isolevel set Tq is the set of matrices X(y) feq. (|2.2[) ) associated with the 
spectral curve C. Now we construct a map from Tq to V\c 9+N ~ 1 (C), called the 
eigenvector mapping, which plays a very important role in the present method. 

Let X = X(y) be an element of Tq. If (x,y) e C, there exists a complex 
iV-vector v(x,y) such that X(y)v(x,y) = xv(x,y), up to constant multiple. 
Then there exists a Zariski open subset C° of C over which the morphism 
C° 9 (x.y) 1 > v(x,y) G P w_1 is uniquely determined. Moreover, for a smooth 
C, this morphism can be extended uniquely over the whole C. Denote this 
morphism by ■ C — > P^ -1 . 

The eigenvector mapping tpc ■ Tc — * Pic d (C) (<i = 5 + A^ — 1) is a map 
defined by the formula: 

W (I) = n(0 P »-i(l)), 

where O p jv-i(1) is the invertible sheaf of hyperplane sections over P^ -1 . Note 
that it is nontrivial to prove ifc (X) e Pic d (C) (see CP §2). 

The role of the eigenvector mapping is to embed the set Tc into Pic d (C). 
The following proposition is originally obtained in van Moerbeke, Mumford 

Proposition 2.1 (|2j, thm. 3) The eigenvector mapping ipe ■ Tc ^ Pic d (C) 
is an embedding. 
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Although the definition of the eigenvector mapping is abstract, we can have 
an explicit formula to express tpc(X) in the present situation. 

9i 

Lemma 2.2 QlJ, §2) Letv(x,y) = | : ] be an eigenvector of X(y) belong- 

ing to x (gi = gi(x 1 y), i = 1, . . . , TV). Then it follows that tpc(X) — [(<7i/<7jv)oo]- 

On the other hand, the divisor (gi/gjv) has the following expression ([2] 
prop. 1): 

(gi/gw) =Vx + {N- l)P -V 2 -(N- 1)Q, (2.4) 

where T>\ and T>% are general and positive divisors of degree g. 

Let Q(X) := V 2 . Lemma I2T21 is rewritten as tp c (X) = [d(X) + (N - 1)Q]. 

2.2 Linearization theorem 

Consider the N x N matrix X t (y) defined by (|2 . 2[) and the associated spectral 
curve C. Let a and r be the isomorphisms on Tq defined by: 

a(X t (y)) = SX t (y)S-\ (i(X t (y)) - J R i (y)X t (y)i? t (y)- 1 - X t+1 (y), (2.5) 
f° 1 \ 

where S = 







For the hpdToda equation Q1.1H1.3| I23 J) , a is the 

1 

\y 0/ 

n-shift operator: n i— * n + 1 and /i is the t-shift operator: t i— > t + 1. 

By calculating the divisors D(cr(X t )) and 5(/i(X t )), we have the following 
theorem which illustrates the flow of the hpdToda equation on Pic d (C): 

Theorem 2.3 (|lj) (I).' Let T> be the divisor T> = P — Q . Then the following 
diagram is commutative. 

r c -> Pic d (c) 

o- I l +[23] . 

T c -» Pic d (C) 

(II); Lef (j = 1,2,..., M) be the divisor £j = P - Aj. If t = j (mod M), 
the following diagram is commutative. 

T c - Pic d (C) 

/i 4 I +&] • 

r c -» Pic d (c) 
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Corollary 2.4 The time evolution 1 1— > t + M is expressed as Z *—> Z + [B — Q] 
onPic d {C). 

Proof. By theorem O (II), on Pic d (C), {!h(| M} is expressed by the 

formula: Zh 2 + [MP — Aq — Ax A M -i]- Then the relation (x) = 

-MP- Q + A Q + At H h + S e Div°(C) yields the result. ■ 



Corollary 2.5 The divisor T)\ in \2.J^\ satisfies T>t = D(a(X t )). 

Proof. By (gU), [V x ] = [d{X t ) + (N - 1)Q — (N — 1)P] - [©(o-^^t))] = 
[t)(cr(X t ))]. Because X>i and t>(cr(X t )) are general, positive and of degree 3, it 
follows that T>i=X>(a(X t )). ■ 
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Corollary 2.6 Let v{x,y) = | : \ be an eigenvector of X(y) belongs to x. 

Then (i) (si/sat) = O(ctA) + (N - 1)P - — (N — 1)Q, and 
(ii) {9n/v9n-x) = i>PO + - 1)^ - - - 

Proof. Part (i) follows immediately from (|2.4j) and corollary |2.5l Applying (i) 
to the matrix a~ 1 X = S^^^XS and noticing that S ■ (gjy V >9i, ■ ■ ■ , gN-i) T — 
(31,32, ■■■ ,9n) T , we obtain (ii). ■ 



Remark 2.1 TTie time evolution 1 1— ► t + M is given by the map: v{X t {y)) := 
LJ 1 {y)X t {y)L t {y). In fact, A2.SH2.3t proves that u(X t (y)) = A t+Af (y). 



3 Tau function solution of the hpdToda equa- 
tion 

In this section, we assume g.c.d.(AT, M) = 1. 
3.1 Construction of tau functions 

We construct a theta function solution of hpdToda equation. As in the previous 
section, X t — X t (y) denotes the square matrix defined by (|2.2p . 

Let C be the (smooth) spectral curve associated with X t . Fix a symplectic 
basis ai, . . . , a g ; flx, . . . , /3 g of C and the normalized holomorphic differentials 
LJi,...,u) g such that J ujj = Sij. The 3x3 matrix f2 := (fg is called 

the period matrix of C. For a fixed point po € C, the Abel-Jacobi mapping 
A : Div(C) — > C S /(Z 9 + f2Z s ) is the homomorphism defined by: 

E y i - E #j ^ E(£ wi, • ■ • , w g ) - xxjj wi, • • • , ug). 
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Let us consider the universal covering tt : il — > C and fix an inclusion i : 
C ^ il. For simplicity, we slightly abuse the notation 'V and "<." to express the 
derived maps Div(il) — > Div(C) and Div(C) «-» Div(il), respectively. Naturally, 
there exists a continuous lift A : Div(il) — > C 9 such that A o t(po) = 0. For the 
projection p : C a -> C a /(Z a + QZ 9 ), it follows that p o A = A o tt. 

For fixed i 6 Z, assume that some lifted positive divisor £>(X t ) e Div 9 (il) 
with 7r(S(Xt)) = fpQ) is specified. Then there uniquely exist two positive 
divisors D(aX t ), T)(pX t ) G Div s (il) such that: 

2(S)((7X t ))=A(a>(Xt) + tP-tQ), Tr(®(aX t ))=U(aX t ), (3.1) 
A(2)(/*X t ))=A(2)(X t )+ t P-tAj) ) 7r(S)(/iXt)) = D(/*X t ), (3.2) 

where t = j (mod M). 

Let r* be a holomorphic function over il defined by the formula: 

T t (p) = 6(A{®(X t )-p-LA}), pell, (3.3) 

where #(•) = 0(»;f2) is the Riemann theta function and A £ div s_ (C) is the 
theta characteristic divisor of C ([3], Chap. II, cor. 3.11). To avoid cumbersome 
notations, we often omit the letters "A", V and use a simpler expression 
r*(p) = 6(D{X t ) — p — A) when there is no confusion possible. 

Although defined over il, r*(p) can also be thought of as a multi- valued 
holomorphic function over C . By the Riemann vanishing theorem ([3], Chap. II, 
thm. 3.11), the zero divisor of r*(p) corresponds with t)(X t ). 

Let r|(p) := 9(D(aX t ) -p — A). Then, by theorem^ the function 

r£(p) ■ r t+1 (p) 0(3>(o* t ) - p - A) • 0(S(/iX t ) — p — A) 



*'(p) := 



r*(p) • r| +1 (p) 0(3) (X 4 ) - p - A) • 6(p(jurX t ) - p - A) 



satisfies [(the zeros of denominator)] = [(the zeros of numerator)] € Pic 2s (C) 
and therefore, it is a single-valued and meromorphic function over C . 



Consider an eigenvector X t (y) 





— X 


(9\\ 






{9%) 



{g\ =9l{x,y) =gf(p)). 



From the relation (gl/gjf) = *{<rX t ) + {N — 1)P - d(X t ) - (N - l)Q (corollary 
2.6|) we derive the following equation by means of Liouville's theorem: 



g\{p)-9T(p) 
gM-gl^ip) 

By virtue of (|3.4I) . we can calculate some special values of 4 r *(p) 



(p) = c x y l KF > %^ , c : constant. (3.4) 



Lemma 3.1 On condition that g.c.d(7V, M) = 1, we ftawe (i) \Ir(P) = c, (ii) 
**(Q)=cx^. 



G 



Proof. The lemma is proved by an elementary calculation, which we shall give 
in the appendix. ■ 

Because 6{D(X) - iQ - A) = 6(5)(X) + (iP - iQ) - iP - A) = 9(D(aX) - 
lP — A), it follows that 

**«) = **+(*)■ where *\ 0») = ^ ' ' 

Then lemma O implies I^+(P) = P N ^\P). 

Repeating this argument for \E' + (p), we also derive ij^Vf- (P) = Il^+(P)i 
and inductively, we have that: 

Let ^ := (n "+"s). Finally we obtain the equations *^ +A r = 

and = where the number d does not depend on n. 

Next consider the following single- valued meromorphic function over C: 

= r t (p)-r t + M (p) = fl(S)(X t )- P -A)-g(S)(^)- P -A) 
W ' r|(p).ri +M (p) e(S(crX t )-p-A)^(S(zy C r- 1 X t )-p-A)" 

Using corollary 12.61 and Liouville's theorem, we derive the following expression: 

<(p)-g^ +M (p) 
9i(p)-9% M i(p)-y' 

which again allows us to compute some special values of $*(p). 



&(p)=c>X y »™ , C: constant, (3.5) 



Lemma 3.2 On condition that g.c.d(iV, M) — 1, we /iaue (i) $ (P) = c', (ii) 



u* 

Proof. See El ■ 

Due to $*(Q) = <S>\(P) and lemma EH we have V£$*_(P) = V^_ 1 $*(P), 
which implies 

^_i$*(P) = v^$V(P) = v^\ + {p) = v$** +++ (p) 

Let := $*++■■■+ ( n "+" s )- Therefore we obtain & n+N = $* t and 

Ki^ti = d', where the number d! does not depend on n. 

Define ri x := r^cP), := r|(tP), r| := r| + ( t P), ■ • • ,ti := r|+...+ (tP) 
(n "+"s). By the arguments above, J* and 1^ have following expressions: 

T t . t+l t . t+M 

4 = dx " , Vt = d>, "+ 1 . (3.6) 
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3.2 Solution of hpdToda 

For ^-dimensional vectors a and b, (a, b) denotes a T b £ C. 

By periodicity d(a N X t ) = d(X t ), there exist integer vectors n, m £ I 9 such 
that A(N(lP — lQ)) = n + Vim. Considering the definition of the Riemann 
theta function (see [3], §11.1, for example), we have 

T n+N = T n x exp(-27ri • (m, z) — 7ri • (m, flm}), i = V—l, 

where z = A(D(a n+1 X t ) — iP — A). By (|3Tgj) . we have 
T t . T *+i 

AA'-^n =d N x -± — ^±^- = d N x exp(-27ri ■ <m,A(tP-tAj)». ( 3 - 7 ) 

T iV+l ' r l 
* r r* • r t+A/ 

k x v 2 v N - a x - t+M 

= d' N x exp(-2i7r • (m, A(lAq + • • • + iA M -i - (M - l)tP- iQ))), (3.8) 

where j = t (mod M). Recall U n 4 +M = II„ 4 and U n = EL V*, which 
imply that d depends on t (mod M) and that d! is independent from t. Finally 
we obtain the conclusion: 

Theorem 3.3 If g.c.d.(JV, M) = 1, I^MKM solves the hpdToda lTOHO|) . 



4 The general cases 

In the previous sections, we have assumed that g.c.d.(7V, M) = 1. Unfortunately, 
the method which we have established in this paper cannot be applied in the 
general cases. 

For example, when N — M — 2, the characteristic polynomial of the matrix 
X t (y) (equation §T2§) is: 

det (X t (y) - xE) = y 2 - y(2x + Ui) + x 2 - U 2 x + U 3 - U A y-\ 

where Ui = / 1 '/1+/J +1 /| +1 +V 1 t y 2 t , U 2 = I t 1 I t 1 +1 +I^I t 2 +1 +I t 1 V 2 t +I t 1 +1 Vl+I t 2 V 1 t + 

i* +1 v$, u 3 = ilP 2 ii +1 it +1 +il +1 it +1 vM+vMnil u 4 = PjUl^il^vfvl 

However, the hungry Toda system (|1 . II 1 1 . 3[) has the extra conserved quantity 
I[ + J| + + 1 2 +1 + VI + V 2 , which is independent from Ui, U 2 , U 3 and [7 4 . 
This means that the spectral curve docs not faithfully reflect the data of the 
system. 

For this reason, we should try to trace the problem to the case g.c.d.(iV, M) = 
1. Denote by Todajv,M the hungry Toda system f| 1 . 1111 . 3 j) associated with the 
positive integers N and M. It is sufficient to prove the following statement. 

Proposition 4.1 Define the initial values 7„ := ( + o(£), (( — > oo,Vn) for 
some complex parameter and let {/*, V^} n- t be a solution o/Todajv,Af- When 
£ — > oo, the new sequence 

rjkM+l jkM+2 rkM+M-1 r/fcM+1 i/fcM+2 v kM+M-l\ 
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is a solution of Toda^M-i ■ 

Proof. We shall prove the following: 



jkM+M-l = jkM-1 + ykM-1 _ ykM+l + ^ (4J) 



jtcM-lykM-1 

K kM+1 = "!L + m-i ■(! + »(!))■ (4-2) 
By <jl.lHl.3p and Remark (page[T]), we have 

P-C + o(C) Nn) => I 7 « +M = C + ° (C) ' (Vn) (C - oo) 

Then, in our situation, it follows that V^ M+1 = V^ M + o(l) for all k e Z> and 
n. Using (|1 .1H1 ,3[) again, we derive equations (|4.1I4.2[) . ■ 



Applying proposition 14.11 repeatedly, we can trace the problem to the case 
g.c.d.(AT,M) = 1. 

Example The hungry Toda system with N — M — 2 can be traced to the case 
N = 2,M = 3. 

Let L := ^ ^ \ R := Jj, i?i := i? 2 := 

1 T i ) ■ Define Xq := LqR2R\Rq. The characteristic polynomial of Xq 



det (X - xE) = -t, 3 + y 2 (C 2 + C/i) - y{(2C + t/ 5 )z + [/iC 2 + C/ 3 } 
+x 2 - (U 2 ( + U 6 )x + U 3 ( 2 + Ui- U A Q 2 y-\ 

where U 5 = 7? + /J? + 1\ + 1\ + V? + V$ and U 6 = I^l\V? + I$I%V£. Note that 
U5 is the hidden conserved quantity o/Toda2,2- Let {In,Vn} n ,t be the solution 
o/Toda2,3. Then the sequence 

Km 7°, lim I*, Km I 3 , Km I*, Km 1%, . . . ; 



lim V° 7 Km V„, lim V^, lim V*, jim V%, . 



/ 6 

£ — >00 £ — >OG C, — >00 C, — >OC £ — >00 



solves Toda2.2- 
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A Proofs of lemmas 



Let \P* (p) and <I>* (p) be the meromorphic functions defined in section [3l We 
shall now prove lemma [37X1 13.21 In the appendix, we assume g.c.d.(N, M) = 1. 

Denote the set of N x TV matrices by A/jy(C) and the subset of diagonal 
matrices by T C M^(C). For a matrix X £ Mjy(C) and subsets A, B C A/jy(C), 
let A+X := {a+X \ a £ A}, AX := {al | a e A}, A+B := {a+b \ a £ A, b £ B} 
and AB := {ab \ a £ A, b £ B}. 

For two meromorphic functions /, 5 over C and a point p £ C, "/ ~ g around 
p" means < lim z ^ p \f(z)/g(z)\ < +oo. 

Let (51,52, ■ ■ • ,9n) T be an eigenvector of X = X(y) £ Tc belonging to an 
eigenvalue x. Then 51, ... , g^ are meromorphic functions over C. The following 
lemma is fundamental. 

Lemma A.l (i) Let k be a local coordinate around P. Then gi/gN = k 1 ^^ 1 + 
■■■ , 92/gN = k N ~ 2 H , . . . , g N -i/g N = k H . 

(ii) Let k be a local coordinate around Q. Then gxjgjss ~ k~ N+1 , g^jgN ~ 
/■• V ' - 9n , 5 V-A- • 

Proof, (i) Recall that we have x = k~ M + ■ ■ ■ and y = k~ N + ■ ■ ■ around P. 
By l|23J) . X t is contained in the subset {E + TS- 1 ){T + S) M = rS- 1 +r + rS' + 
h rS' M_1 + S M . Then the equation X t (y) v = x v implies: 

(-1-iS- 1 + 70 + 71 5 + • • • + 7Af-iS M ~ 1 + S M ) ■ v = k- M v + (higher terms), 

where 7, (i = — 1, 0, . . . , M — 1) are diagonal matrices. Let T := kS. Therefore 

we obtain (t m + ^^f~_i fc M_ *7iT 1 ^ ■ v = v + (higher). Because N and M are 

relatively prime, the solution of Tv = v is v = (fc^ -1 , k N ~ 2 , . . . , 1) T up to a 
constant multiple. This fact leads to the desired result. 

(ii) Let k be a local coordinate around Q such that x — Ek" 1 + ■ ■ ■ and y = 
k M H (Section 12). It follows that 

(7-iS- 1 + 70 + 71 5 + • • • + Jm-iS™- 1 + S M ) ■ v = Ek- X v + (higher). 

Let U := k~ x S. Then we have (j-iU' 1 + J^Zo k l+1 ^U 1 ^ -v = Ev + (higher). 

Standard results from linear algebra prove that there exist (N — 1) complex 
numbers c±, . . . , Cjy-i such that 

U ■ ( Cl k- N+ \ c 2 k~ N+ \ ...,lf = E- ( Cl k- N+ \c 2 k~ N+2 , . . . , If, 

which leads to the desired result. ■ 
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Proof of lemma 13.11 

The equation X t+1 (y)R t (y) = R t (y)X t (y) (EH) implies {g{ + \g\ +1 . . . , g^ 1 ) = 
Rt(y) ■ (gh 92 ■ ■ ■ > 9n)- Tncn gives rise to 



9i I%9n + 9iV 



By lemma EH satisfies = cH , around P, and = c • (7^ + ■ 

around Q. 



Proof of lemma 13.21 

As mentioned in remark [2.11 one has that L t (y)X t+ M{y) — A t (y)L t (y), which 
implies (glgl..., g* N ) = L t (y) ■ (g[ +M \g\ +M . . . , g l + M ). Then leads 



yt t+M , t+M t+M 
v N-l 9n-1 + 9n 9n 



vh9% M y-' + g{ +M gTi-v 



By lemma lA~l~l $* satisfies $' = c'H , around P, and <f>* = c'-(V} j -_ 1 /V^)+- 

around Q. 
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